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A formalism for describing an all-sky map of the polarization of the cosmic microwave background 
is presented. The polarization pattern on the sky can be decomposed into two geometrically distinct 
components. One of these components is not coupled to density inhomogeneities. A nonzero amplitude 
for this component of polarization can only be caused by tensor or vector metric perturbations. This 
allows unambiguous identification of long-wavelength gravity waves or large-scale vortical flows at the 


time of last scattering. 
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With the COBE detection of large-angle anisotropy in 
the cosmic microwave background (CMB), results from 
numerous balloon-borne and ground-based experiments, 
and the advent of a new generation of satellite missions, 
the CMB is becoming an increasingly precise probe of the 
early Universe. CMB anisotropies will help determine 
whether density perturbations (scalar modes) are the 
result of inflation, topological defects, or perhaps some 
other mechanism. Detection of a stochastic gravity- 
wave background (tensor modes) [1] or vortical motions 
in the primeval fluid (vector modes) would help to 
discriminate between these models. Inflation damps out 
vector modes but will produce some tensor modes and 
also predicts a specific relationship between the spectrum 
of the scalar and tensor fluctuations [2]. In contrast, 
topological defects will produce a mixture of scalar, 
vector, and tensor modes. Scalar modes give rise to both 
the observed large-scale structure and CMB fluctuations, 
while in the foreseeable future we can only expect to 
observe the consequences of tensor or vector modes 
through their effects on the CMB. Without a model of 
primordial fluctuations, the contribution of scalar, vector, 
and tensor modes to the CMB temperature anisotropy are 
indistinguishable. 

However, any mechanism which produces temperature 
anisotropies will invariably lead to nonzero polarization 
as well [3-6]. As we demonstrate in this Letter, this 
polarization signal can be used to discriminate between 
scalar and vector or tensor metric perturbations. COBE 
has already mapped the polarization pattern with an 
angular resolution of 7° (although the data have not been 
analyzed), and MAP [7] will measure the polarization 
with a resolution of around 0.3°. 

In prior work, the autocorrelation and cross correlations 
between the Stokes parameters Q and U and the tempera- 
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ture T have been considered. Here Q and U are defined 
with respect to particular orthogonal axes on the celestial 
sphere. While this formalism does provide a complete de- 
scription of the polarization, there is no rotationally invari- 
ant way to lay down orthogonal basis vectors on a sphere, 
so the meaning of QQ, QU, UU, QT, or UT will depend 
on absolute positions of the points being correlated rather 
than just the relative position. Calculations of this sort 
have been done with a small-angle approximation, since 
rotational noninvariance disappears when considering only 
a small patch of the sky. However, this formalism is not 
optimal for describing the complete temperature and po- 
larization correlations present in full-sky maps. 

Here we present a rotationally covariant formalism 
for describing the polarization pattern on a full sky. 
The Stokes parameters, defined by the 2 X 2 correlation 
matrix of the electric field of incoming photons, can be 
described as a tensor field on the celestial sphere. The 
Q and U parameters, describing linear polarization, are 
just given by the symmetric trace-free (STF) part of this 
tensor. For example, in spherical polar coordinates (0, ¢), 
where the spherical metric is gap = diag(1,sin*@), the 
polarization tensor is 


Pâ) = H aa 


—U(ñ)sinð 
—U(f)sind } (u 


— Q (ñ)sin? 0 


The “ab” are the tensor indices, and we use standard 
tensor notation throughout. It is natural to decompose 
the linear-polarization pattern into STF tensor spherical 
harmonics [8,9], which constitute a complete orthonormal 
set of rank-2 STF tensors on the sphere. There are two 
types of harmonic STF tensors, Y, ae and Y ano one of 
each for every one of the usual spherical harmonics Y(jm) 
with Z = 2. Two sets of tensor harmonics are required 
as there are two modes of linear polarization, Q and U. 
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Since Compton scattering can produce no net circular 
polarization, the CMB is expected to have V = 0, and 
the V Stokes parameter will not be considered further. 

The harmonic expansion of an all-sky map of the CMB 
temperature and polarization can be written 


Tâ) 4 
To =] + > alm Yam ô), 
Pa (ii) - @) 
a G G A C E A 
n > laam Y Cna Â) + alim) Yima Â)]. 
lm 
The mode amplitudes are given by 
1 ee 
Aim) = A dn Tâ) Yim Ô), 
1 R 7 ee 
dim) = 7 f dû Pa Â) Yn) Ô), (3) 
1 a X EZA 
any = 7 | di Pa YERA, 


which can be derived from the orthonormality properties 
f dû Yim) (A) Ym (A) = 61 bmn; 


f dn Yeimap AYGA) = Ow Omm' > 


(4) 
f dn Y ras AY Gan Â) _ Ôl O mm! > 


f dû YG ap GY GÂ) = 0. 


Here To is the cosmological mean CMB temperature, and 
we are assuming Q and U are measured in brightness 
temperature units rather than flux units. 

The two geometrically distinct tensor harmonics are 


1 
Y nab = Ni(Y(m):ab = a Bab? Gmieo 3 
N; (5) 
Yim)ab Ea y YamacE + Yum):bc €a) 


Here N; = /2(/ — 2)!/(/ + 2)! is a normalization factor, 
€ap is the completely antisymmetric tensor, and “:” 
indicates a covariant derivative on the sphere. In two 
dimensions, any STF tensor can be uniquely decomposed 
into a part of the form A.g, — (1/2)gapA.¢ and another 
part of the form B.,.€% + B.y-€% where A and B are 
two scalar functions. This decomposition is quite similar 
to the decomposition of a vector field into a part which is 
the gradient of a scalar field and a part which is the curl of 
a vector field; hence we use the notation G for “gradient” 
and C for “curl.” Since the Y(j,)’s provide a complete 
basis for scalar functions on the sphere, the Viejas S and 


Y Emap S provide a complete basis for G-type and C-type 
STF tensors, respectively. This G/C decomposition is 
also known as the scalar/pseudoscalar decomposition [9]. 

In (6, &) coordinates, where Eq. (1) holds, the harmon- 
ics are given explicitly by 


X(im)sind 
—Wom)sin?@ P 
(6) 


= 
Ni f -X Win)Sind 
c aya NI (Im) (Im) 
Ym)ab (fi) 2 ( Wom)sin0 X(imsin? 0 ) 


Ni f W 
G a) l (Im) 
Y(im)ab Â) ( X (im)SinO 


with the definitions 


a a? a. m 

(im) = \ gaz T ClO Be + sino 0 

Xm) = 2in( 2 = cid Ji. ee 
sind \ 30 


The exchange symmetry {Q, U}-—{U,—Q} as GC 
indicates that Yanai and Y(jim)ap represent polarizations 
rotated by 45°. 

A most useful property of the G/C decomposition is 
that, in linear theory, scalar perturbations can produce only 
G-type polarization and not C-type polarization. This is 
in contrast to tensor or vector metric perturbations which 
will produce a mixture of both types. To understand 
why scalar metric perturbations do not produce a C-type 
polarization pattern, consider a scalar perturbation with 
single Fourier mode k in the ĉ direction. The polarization 
in a given direction can be represented by a magnitude 
P = (Q? + U?)!/2 and an orientation angle a from the 
axis defining the Stokes parameters (here, choose 6), where 
tan2a = U/Q. For scalar perturbations, the orientation 
of the polarization can be determined only by the direction 
of k: thus a = 0 if the polarization orientation is along 
the direction of k, or œ = 7/2 if the orientation is 
perpendicular to the direction of k. In either case, in a 
given region of the sky all of the orientations are parallel 
and thus the polarization pattern has no curl. Since the 
curl is a linear operator, summing over Fourier modes 
does not alter this conclusion. For tensor and vector 
perturbations, the azimuthal symmetry in the scalar case 
is explicitly broken, and thus the Fourier vector does 
not completely define the direction of the polarization 
orientation. Another way to state this argument is that 
scalar perturbations have no handedness so they cannot 
produce any curl, whereas vector and tensor perturbations 
do have a handedness and therefore can. 

Finding a nonzero component of C-type polarization in 
the CMB would provide compelling evidence for signifi- 
cant contribution of either vector or tensor perturbations 
at the time of last scattering. Given a polarization map 
of even a small part of the sky one could in principle 
test for vector or tensor contribution by computing the 
combination of derivatives of the polarization given by 
P$? ef, which will be nonzero only for C-type polariza- 
tion, i.e., when vector or tensor perturbations are present. 
Similarly only G-type polarization contributes to P“’.,y. 
Of course, taking derivatives of noisy data is problematic; 
more robust measures are given below. 

We now turn to statistics of CMB polarization. If the 
cosmological inhomogeneities are Gaussian random noise, 
then to the extent linear theory is valid, the CMB fluctua- 
tions will also be Gaussian random noise. Regardless of 
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whether the distribution is Gaussian, rotational invariance 
requires that the two-point correlations be of the form 


(alin) Alm’) a CT 8 Sim! » 
(atin) Atm’) = CP SW Simm’ » 
(atin) tim’) = Cf ôw b mm’, 
(aim) Atm’) x CTS Simm’ > 
(alin) Atm) E cs 81 Ômm' > 
(atm agim) z Cr 61 Ômm . 


If we also require that the distribution of inhomogeneities 
be invariant under parity, then Ge = cre =0 since the 
Yaim) and the ‘ea have parity (—1)! while the Yonint 
yet, 


(8) 


have parity (—1 Measuring a nonzero Cr and/or 
CPF would be quite interesting, indicating a handedness 
to the inhomogeneities in our universe. However, we 
do not expect this and will henceforth only consider the 
four angular power spectra {C/, C PC, CPS}. The first 
is the well-known angular power spectrum of temperature 
anisotropies while the last three, new to this paper, are re- 
lated to various quantities in previous work (see Ref. [10]). 
Note that the scalar, vector, and tensor contribution to each 
of the C;’s adds in quadrature, i.e., for X = T, G, C, TG 


(are — cot Ẹ Cee he Cee. (9) 
and this is true whether or not the fluctuations are Gaussian. 
We have argued that CE" = 0. 


Given an all-sky temperature-polarization map, one can 
determine the aqm) s using Eq. (3), and then construct 
estimators for the C;’s in the usual way, i.e., 


~ L laq jl? ~ l lati iI? 
Gey Sel Ge yy Fer 
KEEN 214+ 1 Meran I+1 
‘oS L lak? a i „Ts G (10) 
cE = Am) cis = 5 Gm) 4(Im) 
Re) ee i i Sı 2+1 


If only part of the sky is mapped, the same techniques 
developed to analyze anisotropy with incomplete sky 
coverage [11] may be applied to polarization to construct 
other estimators of the various C;’s. The mean square 
polarization is 


Q? + U2=2P%P,, = PE + PÈ, (11) 


where 
Pe = U+14 Pe +14 
= Cr, = Cr. 
To 2 87 pe 2 8T 
(12) 


Since scalar modes do not contribute to P , this statistic 
provides a powerful and unambiguous model-independent 
probe of tensor and vector perturbations. 

To test a given spectrum of tensor modes against a po- 
larization map, comparing the complete set of predicted 


c$ with the estimators cr is more powerful than con- 
sidering only PŽ, if the detection has sufficient signal 
to noise. Usually, however, the theory being tested has 


2060 


scalar as well as nonscalar modes, along with undeter- 
mined cosmological parameters. If so, the most infor- 
mation can be extracted from the map by comparing the 
entire set of predicted moments, {C E CS, ang CP. with 
the measured estimators [12,13]. 

Note that only the C rs potentially allow detection of a 
small vector or tensor signal. If scalar perturbations domi- 
nate, then the vector or tensor signal in {C7}, C$, c79} 
may be swamped by the cosmic variance in the scalar sig- 
nal, but the Cy ’s are not contaminated in this way. The 
cross-correlation moments Cre, which differ for scalar, 
vector, and tensor perturbations [14], will be larger than 
the polarization autocorrelation moments. Therefore, the 
temperature-polarization cross correlation may be mea- 
sured with some precision. 

Much of the small-angle formalism of Refs. [3-5,14- 
17] can be reproduced by replacing the Yam) (ñ)’s in our 
formalism with Fourier modes, e", and using regular 
derivatives rather than covariant ones. This small-angle 
formalism is completely analogous to that developed 
above and will provide an accurate description of a region 
of sky small enough to be approximated by a flat surface. 
The G/C decomposition in the small-angle formalism can 
be used to detect nonscalar perturbations on small scales, 
though the tensor and vector signal are liable to drop off 
rapidly at angular scales smaller than a few degrees. 

To make contact with previous work, we can write the 
two-point temperature and polarization correlation func- 
tions [3-5,14-17] in terms of multipole moments [10]. 
Although correlation functions of Stokes parameters which 
appear in the previous literature depend on the positions 
of the points being correlated, rotationally invariant cor- 
relation functions exist which are closely related to those 
discussed above. To construct them, define Stokes param- 
eters Q, and U, with respect to axes which are parallel and 
perpendicular to the great arc (or geodesic) which connects 
the two points being correlated. The two-point correlation 
functions are 


C" (0) 


To (T ÂT (iz) A, =cos 8 


21+ 1 
= 5 i ct P}(cos@), 
; T 


T) (Q, ûı)Q, (fi2))a, *fin=cos 0 


21+ 1 
5y 5 MLC? Waz (0,0) 
7 T 


— iCF X(2)(0, 0)], 
To (U, (fi) U,-(fi2)) A, -fts—cos 6 


21+ 1 
2 5 Ni[CF Wan (0,0) 


I T 


c2 (0) 


c” (0) 


— iCF Xa (0,0)], 
C72 (0) = Tg” (T (1) O, (fi2)) A, cos 6 


+1 
=> Ni CTS P?(cosd), (13) 
T 4m 
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where @ is the angle separating the two points and P? is the 
m = 2 associated Legendre function. Since Q, and T are 
invariant under reflection along the great arc connecting the 
two points while U, changes sign, (0,U,) = (U,T) = 0 
if statistical invariance under parity holds. Equations (13) 
reduce to the correct small-angle formulas [17] when 0 « 
1. The functions {C7(6), C?(0), U (0), CT2(0)} are a 
different way of representing {C}, C7, Cr, CFO} and vice 
versa. In Gaussian models either set provides a complete 
statistical description of the temperature and polarization 
patterns. 

The largest hurdles to detecting and characterizing 
CMB polarization are sensitivity and foregrounds. In 
adiabatic models with standard recombination the po- 
larization is only a few percent of the anisotropy, al- 
though it may be larger in reionized [16], isocurvature, 
or topological-defect models. Thus the polarization sig- 
nal is at least an order of magnitude below current 
experimental sensitivities. Experiments planned or en- 
visioned over the coming decade, however, will likely 
attain the raw sensitivity necessary for detailed polariza- 
tion investigations. Polarized emission from foreground 
sources is a relatively unknown factor at this time. Fore- 
ground emission and any Faraday rotation will certainly 
contribute to the C-type polarization, but these contami- 
nants can be subtracted using multifrequency observa- 
tions. On subdegree scales, where the signal from vector 
and tensor modes are liable to be negligible, any mea- 
surable C-type polarization is a likely indicator of con- 
tamination. Polarization measurements will be difficult, 
but the promise of using them to detect gravity waves or 
vorticity, and hence to discriminate between cosmologi- 
cal models, makes these measurements potentially very 
valuable. 
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